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Abstract 

We study the influence of ferromagnetic and antiferromagnetic bond defects on the ground-state energy of antiferromagnetic spin 
chains. In the absence of translational invariance, the energy spectrum of the full Hamiltonian is obtained numerically, by an iterative 
modification of the power algorithm. In parallel, approximate analytical energies are obtained from a local-bond approximation, 
proposed here. This approximation results in significant improvement upon the mean-field approximation, at negligible extra 
computational effort. 
© 2008 Elsevier B.V. All rights reserved. 
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1. Introduction 

Many different types of magnetic systems are commonly 
modeled by spin Hamiltonians, such as the Heisenberg 
model, fn special limits, such model Hamiltonians can be 
solved analytically, e.g., by the Bethe Ansatz for integrable 
one-dimensional systems. Away from special limits, power- 
ful numerical methods, such as Monte Carlo simulations or 
density-matrix renormalization group techniques, provide 
much useful information, but at high computational cost. 
Full exact diagonalization is even more expensive, and thus 
limited to rather small systems. Quite generally, analytical 
and numerical techniques work best for models in which 
all sites are equivalent, so that translational invariance can 
be used to reduce the complexity of the problem. 

In magnetic nanostructures[l] and molecular magnets [2] , 
the generic many-body features of spin Hamiltonians coex- 
ist with additional real- life complications, such as impuri- 
ties, defects, boundaries, textures, etc. The resulting spa- 
tially inhomogeneous spin Hamiltonians do not have trans- 
lational symmetry, and often also lack other simplifying 
symmetries, which makes them hard to treat by traditional 
analytical and numerical methods. 

Mean-field theory can, of course, always be applied, but 
is not sufficiently reliable to permit quantitative modeling, 
and is often even qualitatively wrong. In this paper we de- 
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scribe one analytical and one numerical method to obtain 
beyond-mean-field energies for spin chains without trans- 
lational invariance. 

The analytical calculations are based on model density- 
functional theory (DFT). In ab initio electronic-structure 
calculations DFT [3] is a useful way to include correlations 
beyond the mean-field approximations, at very little addi- 
tional computational cost. To deal with defects within DFT 
for spin chains, we propose a local-bond approximation 
(LB A), akin to the local-density approximation (LDA) of 
ab initio DFT and the local-spin approximation (LSA) pre- 
viously proposed for spin systems with impurities[4,5,6,7]. 
These schemes are described in Sec. 2. 

To obtain results with higher precision than is possible 
by local approximations, and to judge the performance of 
such simple improvements on mean-field theory, we also 
perform fully numerical calculations, employing an itera- 
tive modification of the power algorithm, described in Sec- 
tion 3. Representative results are discussed in Section 4. 

2. Local approximations for spin Hamiltonians 

DFT has been applied to model Hamiltonians, such as 
the Heisenberg[4,5,6,7] and the Hubbard[8] model, within 
generalized local-density approximations. Specifically for 
the Heisenberg model with impurity spins, this scheme con- 
sists in adding to the mean-field energy a local approxima- 
tion to the correlation energy E c , of the form 
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(i) 



where the sum runs over all sites i, and e^ om (J, S) is the per- 
site correlation energy of the homogeneous spin S Heisen- 
berg model, for which various approximations exist [4,5]. 
This so-called local-spin approximation [4,5] has been ap- 
plied to investigate the energetics of impurities in one, two 
and three-dimensional Heisenberg models[4,5,6,7]. Impuri- 
ties here are defined, relative to a homogeneous model in 
which all sites have the same spin S, as sites with a spin 
Si ^ S (see inset of Fig. 1 for an example). 

Another important class of spatial inhomogeneity, de- 
fects, cannot be treated within the LSA. We define defects, 
again with respect to a homogeneous model, as bonds hav- 
ing strength J D ^ J (see insets of Figs. 1 and 2 for exam- 
ples). To extend the local-approximation scheme to spin 
Hamiltonians with defects, we here propose the local-bond 
approximation (LB A), 



E. 



LBA 



[J b ,S] 



(J,S)\ 



J—>Jb ' 



(2) 



where the sum runs over all bonds b, and for e^ om (J,S) 
we use the same expressions employed in the LSA. Con- 
ceptually, the LBA decomposes the system in bonds, act- 
ing between sites, whereas the LSA decomposes it in sites, 
connected by bonds. Both approximation schemes become 
exact for infinite homogeneous systems, and both can be 
used in analytical calculations, as all one has to do to ob- 
tain corrections to the mean-field energies is to evaluate 
Eqs. (1) or (2) site by site or bond by bond. 

3. Numerical Ground- State Energy 

To obtain energies of higher quality than is possible with 
analytical calculations employing local approximations, we 
resort to a numerical scheme. Even in the presence of impu- 
rity spins and/or bond defects, the Heisenberg Hamiltonian 



H — JjiSj • Si 



+i ? 



(3) 



where Ji is the exchange integral between nearest-neighbor 
spin- vectors S* and Si+i, conserves the z-component of the 
total spin, S z . The Lieb-Mattis theorem guaranties that 
the GS pertains to the subspace of minimum \S Z \. As ba- 
sis vectors we therefore use the set of quantum numbers 
{|toi, m2, tun)}, where rrii are the eigenvalues of Si, and 
can take values —Si,—Si + l,...,Si. A practical way of gen- 
erating this set of states is by decomposing integer numbers. 
For instance, for a chain with four spins 1 /2, the binary de- 
composition of the integer 6 gives the sequence 0110, rep- 
resenting the vector | — 1/2,1/2,1/2,-1/2) , one among 
other five of the subspace with S z — 0. It is straightforward 
to assemble the matrix representation of H in this base. 
However, although the resulting matrix is sparse, storing 
it is impractical: for 20 spins 1/2 the order of the matrix is 
184756, with 1108536 non- vanishing elements. In practice, 



we therefore do not store the matrix, but compute each 
matrix element every time it is needed; this saves a lot of 
memory, but is time consuming. 

To obtain the ground-state energy of H we propose an 
iterative modification of the well-known power algorithm. 
This modification, to be described in more detail in a sep- 
arate publication, is easy to code and requires less memory 
than the Lanczos method, although it typically takes more 
processing time to extract the ground state. The power 
algorithm starts by decomposing a trial function for the 
ground-state eigenvector in terms of the unknown eigen- 
vectors \ip n ) of H, according to Vt) = X)n a «IV'n), with 
a n constants. If e is an upper limit of the energy spectrum, 
acting k times with the operator H — e on |V>t) yields 



{H-e) k \i) T )=Y j an{E n -e) k \4> n ) 



(4) 



where H\tp n ) = E n \tp n ). For antiferromagnetic chains, the 
highest energy e corresponds to the ferromagnetic configu- 
ration, whose value is trivial even in the presence of impu- 
rities or defects. For k — > oo, the above series is dominated 
by the ground-state term ao(E — e) k \ipo)- Therefore, we 
can extract E by performing 



Eq 



lim 

k— >oo 



{<j>\{H - e) k \i> T ) 



(5) 



where \<f>) is any vector non-orthogonal to the ground state. 
Without the constant e, the method yields the highest (in 
modulus) eigenvalue instead. Our experience shows that 
the best trial function is the Neel state, a mean-field ap- 
proximation for the ground state, consisting of a sequence 
of an up and down spins. 

To speed up the search for Eq, we implement the above 
limit iteratively: at each step k , we use for \<f>) the state 
\<j>k-i) — (ff~ e ) fe_1 IV't) obtained in the previous step. This 
reduces the processing time to reach convergence, which 
we characterize by two successive values of Eq/N differing 
by less than 10~ 13 . Typically, a few hundred /c-steps are 
required for chains larger than 20 spins, and only a few 
dozens for smaller chains. Using a desktop microcomputer 
with 1.5 MB of RAM, we obtained the ground-state energy 
of a homogeneous chain with 30 spins 1/2, reproducing the 
results of Ref . [9] . Bond defects do not require more memory 
space. 

4. Antiferromagnetic spin chains with ferro- and 
antiferromagnetic defects 

As a first application of the LBA concept we compare, 
in Fig. 1, an antiferromagnetic (AFM) spin 1/2 ring (peri- 
odic boundary conditions) with one impurity spin Si = 3/2 
to an AFM spin 1/2 ring with one bond defect Jd = 5 J. 
In the mean-field approximation, both systems are, erro- 
neously, predicted to have the same ground-state energy. 
This spurious degeneracy is lifted by adding the LSA and 
LBA correlation energies, respectively. Judged by the re- 
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maining distance to the exact data, LSA performs slightly 
better for the impurity than LBA does for the defect, but 
both provide significant improvements on the mean-field 
data. 
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Fig. 1. Ground-state energy of an AFM spin 1/2 ring with a spin 3/2 
impurity (triangles), and with an AFM defect of strength Jo = 5 J 
(circles), treated exactly (filled symbols) and via local approxima- 
tions LSA/LBA (open symbols). Within mean-field theory (stars), 
the impurity and the defect yield the same ground-state energy. In- 
clusion of correlation energy in the LSA (impurity case, open trian- 
gles) and the LBA (defect case, open circles) removes this spurious 
degeneracy. 

Ferromagnetic (FM) defects, with Jo < 0, can be han- 
dled similarly. However, an ambiguity arises in how the lo- 
cal substitution is to be performed. In principle, the LBA 
can be applied to a FM defect by substituting J — > \Ji\, 
J — > Jj or J — > 5 j. I j. | . The first possibility can lead to 
positive correlation energies, in violation of the variational 
principle. The second and third possibilities correctly pre- 
dict negative correlation energies, but numerically the sec- 
ond is slightly inferior to the third for the type of system in- 
vestigated here, as judged by comparison to exact data. Be- 
low, we thus employ the third substitution. Physically, this 
choice, J — > means that a ferromagnetic bond does 

not contribute to the correlation energy, which is rather 
reasonable, as at zero temperature the mean-field energy 
of a homogeneous ferromagnetic chain is already exact. 

In Fig. 2 we show ground-state energies of an AFM spin 
chain (open boundary conditions) with one AFM defect 
Jd = +3J, and of the same chain with one relaxed FM de- 
fect, Jd = — 3 J. Comparison of LBA data with exact data 
shows that the LBA significantly improves on the mean- 
field approximation both for FM and AFM defects. Quan- 
titatively, the performance for the FM defect (squares in 
Fig. 2) is better than that for the AFM defect, but in both 
cases a significant improvement over the mean-field data 
is achieved. Note, in particular, that for an FM defect the 
mean-field curve is even qualitatively wrong, predicting a 
wrong sign for the slope at N —* 0, whereas the LBA re- 
covers the correct behavior. 
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Fig. 2. Ground-state energy of an AFM spin 1/2 chain with an 
AFM defect of strength Jr> = +3 J (circles), and with an FM defect 
of strength Jr> = — 3J (squares). Exact data (full symbols) are 
compared to LBA data (open symbols) and mean-field data (stars). 

We stress that the LBA calculations for both AFM and 
FM defects, as well as LSA calculations for impurities, can 
be done analytically, in any dimensionality and for any sys- 
tem size, regardless of boundary conditions and symme- 
tries. Such local approximations thus provide a very con- 
venient way to obtain beyond-mean-field results at the ex- 
pense of a conventional mean-field calculation. When used 
in this way, local approximations yield robust results of 
moderate precision, even for large and complex systems. 

High precision results, on the other hand, cannot be ex- 
pected from simple local approximations. Improved algo- 
rithms for fully numerical diagonalization, such as that de- 
scribed in Section 3, can be used to make progress from 
exactly the opposite starting point: high precision for small 
systems. 
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